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Introduction 

Lysophosphatidic  acid  (LPA)  elicits  a  variety  of  responses  including  mitogenesis,!'^ 
cytoskeletal  changes, 4”  10  activation  of  Ca2+  transients,!  1-13  effects  on  apoptosis.!^  These 

responses  are  elicited  via  the  LPAi/EDG2,  LPA2/EDG4,  and  LPA3/EDG7  G  protein-coupled 
receptors.  These  receptors  are  members  of  the  endothelial  differentiation  gene  family.  In  order 
to  understand  the  physiological  significance  of  LPA  highly  selective  antagonists  are  necessary. 
The  work  of  Imamura  and  colleagues  has  established  that  cancer  cells  require  serum 

factors  for  invasion^  and  that  LPA  is  the  major  serum  factor  responsible  for  this  effect.  This 
invasive  effect  is  mediated  by  the  Rho  signaling  pathway.  LPA  is  fully  capable  of  restoring 
tumor  cell  invasion  in  serum  free  systems  10,16^  The  central  goal  of  the  present  study  is  to 
develop  and  validate  an  LPA3/EDG7  model  and  based  on  this  model  rationally  design  ligands 
with  receptor  subtype-specific  action,  that  will  provide  critically  needed  new  tools  and 
knowledge  in  this  developing  field.  The  development  of  validated  computational  receptor 
models  will  aid  future  attempts  to  develop  such  compounds,  which  may  be  useful  in 
understanding  and  eventually  treating  prostate  cancer. 

BODY 

Three  antagonists  selective  for  LPAi  and  LPA3  have  recently  been  identified.  First, 
DGPP  (8:0)  (Figure  1)  was  shown  to  be  a  selective  inhibitor  of  LPAi  and  LPA3.  DGPP 
inhibited  LPA-induced  Ca2+  transients  in  stably-transfected  RH7777  cells  expressing  LPAi  and 

LPA3  with  K;  of  6600+680  and  106+28  nM  respectively.  17  However  DGPP  had  no  effect  on 
LPA-induced  Ca2+  response  in  stably-transfected  RH7777  cells  expressing  LPA2.  Second 
compound,  FAP  (12:0)  (Figure  1)  was  found  to  be  a  selective  antagonist  of  LPAi  and  LPA3. 
FAP  inhibited  LPA-induced  Ca2+  responses  in  RH7777  cell  expressing  LPA]  and  LPA3  with  a  Ki 
values  of  9  pM  and  89  nM  respectively.  FAP  did  not  inhibit  LPA-induced  Ca2+  mobilization  in 
RH7777  cells  expressing  LPA2,  however  FAP  activated  Ca2+  mobilization  in  these  cells  with  an 
EC50  700  nM  18.  Finally  NAEPA  (Figure  1)  was  also  shown  to  be  a 

selective  antagonist  of  LPAi  and  LPA3.  NAEPA  inhibited  LPA  induced  GTP  [y35S]  binding  in 
RH7777  cells  expressing  LPAi  and  LPA3  with  Ki  of  1 3  7  nM  and  428  nM  respectively.  1 9 
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1-C16-GP(R) 


DGPP  8:0  FAP  12:0 


NAEPA 

Figurel:  Agonist  and  antagonist  of  LPA  receptors 

Inactive  and  Active  models  of  the  LPA  receptors  were  developed  in  our  research  group 
20  An  active  model  of  each  LPA  receptor  was  developed  via  homology  modeling  in  MOE 


5 


program  and  was  based  on  the  validated  model  of  SIPi  22.  An  inactive  model  of  each  LPA 
receptors  was  developed  by  homology  modeling  in  the  MOE  program  and  was  based  on  the 
bovine  rhodopsin  crystal  structure.2^  Autodock  3.02^525  was  used  to  assess  the  docked  energies 
of  DGPP  (8:0)  and  FAP  (12:0),  with  the  inactive  models  of  LPAi/EDG2,  LPA2/EDG4,  and 
LPA3/EDG7.  Both  the  R  and  S  isomers  of  LPA  (18:1)  were  docked  against  active  model  of 
LPA3/EDG7  to  assess  the  stereochemical  preference  towards  the  receptor. 

We  explored  the  binding  site  and  assessed  the  critical  amino  acid  interactions  of 
antagonists  and  agonists  by  docking  to  the  inactive  and  active  form  of  the  receptor  respectively. 
The  antagonists,  DGPP  (8:0)  and  FAP  (12:0),  were  docked  against  the  inactive  receptor  models. 
The  agonists,  R-LPA  (18:1)  and  S-LPA  (18:1)  were  docked  against  the  active  receptor  models. 
One  hundred  complexes  of  receptor-ligand  pair  were  generated,  in  order  to  evaluate  the  binding 
region  of  the  ligand.  The  best  ligand-protein  complex  for  each  ligand  was  energy  minimized 
with  the  MMFF94  forcefield26  to  a  RMSG  of  0.01  kcal/mol-  A.  The  ligands  were  then  removed 
from  the  minimized  receptor-ligand  complexes  and  100  additional  complexes  were  generated 
with  Autodock  in  order  to  re-examine  the  binding  energy  after  allowing  the  receptor  to 
acclimatize  to  the  ligand.  The  docking  runs  were  repeated  until  the  final  docked  energy 
converged. 

We  have  docked  DGPP  (8:0)  and  FAP  (12:0)  in  our  inactive  LPAi,  LPA2,  and  LPA3 
receptor  models  and  our  docked  energies  agree  with  the  observed  trend  in  inhibition  constants 
(Kj).  The  docked  antagonist  positions  overlap  the  agonist  position  only  at  the  polar  head  group, 
whereas  the  hydrophobic  tails  diverge.  In  the  case  of  the  agonist  the  hydrophobic  tail  is  in  the 
transmembrane  (TM)  region  where  as  in  the  antagonist  the  hydrophobic  tail(s)  are  in  the  region 
just  outside  the  helices.  DGPP  (8:0)  occupies  a  position  just  outside  the  helical  bundle  towards 
the  extracellular  end.  The  polar  head  group  of  DGPP  (8:0)  forms  ion-pairs  with  K95 
(extracellular  loop  1),  R105  (TM3),  and  R276  (TM7).  The  polar  head  group  is  also  in  a  position 
to  form  a  hydrogen  bond  with  N172  (extracellular  loop  2),  whereas  the  hydrophobic  tail  is 
oriented  towards  helices  IV  and  V  (Figure  2).  FAP  (12:0)  also  occupies  an  area  outside  the 
helical  bundle  towards  the  extracellular  end.  The  phosphate  group  forms  ion-pairs  with  residues 
R105  (TM3),  R276  (TM7),  and  K95  (extracellular  loopl).  The  hydrophobic  chain  is  oriented 
towards  helices  IV  and  V  (Figure  3) 
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Figure  2:  On  the  left  is  a  ribbon  model  of  the  inactive  form  of  LPA3/EDG7  receptor  with  DGPP 
(8:0).  On  the  right  is  the  close-up  view  of  DGPP  (8:0)  interactions  with  amino  acid  residues  of 
LPA3/EDG7 


The  R  and  S  isomers  of  LPA  (18:1)  were  docked  against  the  active  model  of  LPA3  to 
assess  the  stereo-selectivity  of  the  receptor.  A  previous  study  showed  that  the  R  and  S  isomers 
of  1-C16-GP  (Figure  1)  were  equally  active  in  mobilizing  Ca2+  in  RH7777  cells  stably 
expressing  LPAi,  LPA2,  and  LPA3.27  Based  on  our  results  as  shown  in  Table  1  we  do  not  see 
any  differences  in  energy  with  respect  to  R  and  S  isomer.  On  assessing  the  complex  of  the  R 
isomer  of  LPA  (18:1)  with  the  active  form  of  LPA3  receptor  we  observe  that  the  polar  phosphate 
functionality  is  in  close  proximity  to  ion-pair  with  R105  and  R276.  We  also  observe  that  the 
hydroxyl  group  at  the  sn-2  position  is  in  a  position  to  form  hydrogen  bond  with  Q106  as  shown 
in  figure  4. 
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Figure  3:  On  the  left  is  a  ribbon  model  of  the  inactive  form  of  LPA3/EDG7  receptor  with  FAP  (12:0). 
On  the  right  is  the  close-up  view  of  FAP  (12:0)  interactions  with  amino  acid  residues  of  LPA3/EDG7 


LPA  (18:1) 

Final  Docked 
Energy 
(kcal/mol) 

Intermolecular 
energy  (kcal/mol) 

Torsional  Free 
energy 
(kcal/mol) 

Internal  energy 
of  the  ligand 
(kcal/mol) 

R-isomer 

-15.67 

-16.64 

+6.85 

+0.97 

S-isomer 

-14.78 

-18.13 

+6.85 

+3.35 

Table  1:  Docked  energies  of  R  and  S  isomers  of  LPA  (18:1)  against  LPA3/EDG7 


Figure  4:  On  the  left  is  a  ribbon  model  of  the  active  form  of  LPA3/EDG7  receptor  with  LPA  (18:1).  On 
the  right  is  the  close-up  view  of  LPA  (18:1)  interactions  with  amino  acid  residues  of  LPA3/EDG7 


Key  Research  Accomplishments 

•  Refined  models  of  inactive  forms  of  LPAi,  LPA2,  and  LPA3  G-protein  coupled  receptors 

•  Refined  models  of  active  forms  of  LPAi,  LPA2,  and  LPA3  G-protein  coupled  receptors 

•  Identification  of  a  probable  site  for  binding  of  LPA  inhibitors  DGPP  (8:0)  and  FAP 
(12:0)  to  LPAi,  LPA2,  and  LPA3  using  Autodock  3.0 

•  Identification  of  critical  amino  acids  residues  that  the  above  mentioned  inhibitors  interact 
with 

•  Docked  R  and  S  isomers  of  LPA  (18:1)  to  the  active  form  of  LPA3  receptor  and  our  data 
suggests  that  there  is  no  significant  difference  between  the  binding  of  the  isomers  to  the 
receptor 
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Reportable  outcomes 

•  Vineet  M.  Sardar,  Debra  L.  Bautista,  David  J.  Fischer,  Kazuaki  Yokoyama,  Nora 
Nusser,  Tamas  Virag,  De-an  Wang,  Daniel  L.  Baker,  Gabor  Tigyi,  and  Abby  L.  Parrill. 
Molecular  Basis  for  Lysophosphatidic  Acid  Receptor  Antagonist  Selectivity.  Biochimica 
et  Biophysica  Acta  1582  (2002)  309-317  (Paper  attached  in  Appendices) 

•  Tamas  Virag,  Don  B.  Elrod,  Karoly  Liliom,  Vineet  M.  Sardar,  Abby  L.  Parrill,  Kazuaki 
Yokoyama,  Gangadhar  Durgam,  Wenlin  Deng,  Duane  D.  Miller,  and  Gabor  Tigyi.  Fatty 
Alcohol  Phosphates  are  Subtype-Selective  Agonist  and  Antagonists  of  LPA  Receptors. 

( Molecular  Pharmacology  63(5),  (2003)  in  press) 

•  Vineet  M.  Sardar,  Tamas  Virag,  David  J.  Fischer,  Don  Elrod,  Debra  L.  Bautista,  De-an 
Wang,  Nora  Nusser,  Kazuaki  Yokoyama,  Daniel  L.  Baker,  Duane  D.  Miller,  Gabor 
Tigyi,  and  Abby  L.  Parrill.  Molecular  Modeling  of  Lysophosphatidic  Acid  Receptor 
Antagonists.  ACS  224th  National  Meeting,  Boston,  MA,  August  18-22, 2002  (Abstract 
attached  in  Appendices  MEDI  #  79) 

Conclusions 

We  have  refined  active  and  inactive  models  of  LPAi/EDG2,  LPA2/EDG4,  and 
LPA3/EDG7  receptors.  A  probable  binding  site  just  outside  the  helical  domain  for  LPA3/EDG7 
inhibitors  has  been  identified.  The  amino  acids  that  interact  with  receptor  specific  inhibitors, 
DGPP  (8:0)  and  FAP  (12:0)  have  provided  us  with  useful  insight.  This  information  can  be  used 
to  design  more  potent  and  selective  inhibitors  for  the  prostate-expressed  LPA3/EDG7  receptor. 
Rationally  designing  ligands  based  on  the  ligand-protein  interactions  will  aid  future  attempts  to 
develop  receptor-specific  compounds,  which  may  be  useful  in  understanding  and  eventually 
treating  prostate  cancer  and  restricting  its  invasiveness. 
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Abstract 

Recent  characterization  of  lysophosphatidic  acid  (LPA)  receptors  has  made  possible  studies  elucidating  the  structure- activity 
relationships  (SAR)  for  agonist  activity  at  individual  receptors.  Additionally,  the  availability  of  these  receptors  has  allowed  the  identification 
of  antagonists  of  LPA-induced  effects.  Two  receptor-subtype  selective  LPA  receptor  antagonists,  one  selective  for  the  LPA,/EDG2  receptor  (a 
benzyl-4-oxybenzyl  JV-acyl  ethanolamide  phosphate,  NAEPA,  derivative)  and  the  other  selective  for  the  LPA3/EDG7  receptor  (diacylglycerol 
pyrophosphate,  DGPP,  8:0),  have  recently  been  reported.  The  receptor  SAR  for  both  agonists  and  antagonists  are  reviewed,  and  the  molecular 
basis  for  the  difference  between  agonism  and  antagonism  as  well  as  for  receptor-subtype  antagonist  selectivity  identified  by  molecular 
modeling  is  described.  The  implications  of  the  newly  available  receptor-subtype  selective  antagonists  are  also  discussed.  ©  2002  Elsevier 
Science  B.V.  All  rights  reserved. 

Keywords:  G  protein-coupled  receptor.  Molecular  modeling;  Endothelial  differentiation  gene  receptor;  Diacylglycerol  pyrophosphate;  Homology  modeling; 
Lysophosphatidic  acid 


1.  Introduction 

During  the  last  decade,  enormous  progress  has  been 
made  in  identifying  and  characterizing  the  cellular  receptors 
specific  for  the  phospholipid  growth  factor  lysophosphatidic 
acid  (LPA)  [1-6].  Three  of  these  receptors,  LPA!/EDG2, 
LPA2/EDG4,  and  LPA3/EDG7,  are  members  of  the  endo¬ 
thelial  differentiation  gene  (EDG)  family  and  share  30-35% 
homology  with  sphingosine-1 -phosphate  (SIP)  receptors  in 
the  same  family  (SlPj/EDGl,  S1P2/EDG5,  S1P3/EDG3, 
SIP4/EDG6,  SIP5/EDG8).  Research  into  the  roles  of  each 
of  these  receptors  in  the  biological  functions  of  LPA  has 
consequently  developed  with  groups  examining  knockout 
mouse  models  [7],  expression  patterns,  and  functional  roles 
in  overexpression  systems  (for  reviews,  see  Refs.  [8-11]). 
Until  very  recently,  a  suite  of  receptor-selective  ligands, 
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both  agonists  and  antagonists  that  selectively  activate  or 
block  specific  receptors  has  been  lacking.  Here  we  describe 
progress  that  has  been  made  in  elucidating  the  structure - 
activity  relationships  (SAR)  of  these  three  receptors,  as  well 
as  in  die  development  of  receptor-specific  antagonists. 


2.  LPA  receptor  agonist  SAR 

Early  studies  done  to  elucidate  the  structural  features  of 
LPA  that  determine  its  biological  activity  were  performed  in 
many  cell  lines,  as  the  receptors  responsible  for  the  bio¬ 
logical  effects  of  LPA  were  not  yet  known.  Only  after  the  G 
protein-coupled  receptors  involved  in  LPA  signaling  were 
cloned  and  characterized  [1,2, 4-6]  were  SAR  studies  per¬ 
formed  in  cell  lines  overexpressing  individual  receptors. 
These  data  can  now  be  pooled  and  compared  to  draw 
conclusions  about  the  receptors  responsible  for  the  bio¬ 
logical  end  points  observed  in  various  cell  lines. 

The  modifications  made  to  LPA  can  be  discussed  based 
on  its  structural  regions  identified  as  X,  R1,  R2  and  Q  in  Fig. 
1 .  Of  these  regions,  modifications  to  X  seem  to  completely 
eliminate  activity.  Example  modifications  that  significantly 
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tiosphatidic  add  (LPA  18:0) 

Lysop 

o. 
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Cydic-phosphatidfc  acid  (cPA) 


o 


Diacylglycerolpyrophosphate  (DGPP  8:0) 


Benzyl-4-oxyb  enzyl  N-acyl  etbanolamide  phosphate  (NAEPA) 


Fig.  1 .  LPA  and  analogues  utilized  to  describe  the  LPA  receptor  SAR. 


alter  the  size  of  the  headgroup  include  addition  of  inositol, 
ethanolamine,  choline  and  serine,  which  were  all  inactive 
against  LPAi  or  LPA2-expressing  Tag-Jurkat  cells  [12], 
LPA3-expressing  HEK293T  cells  [2]  or  LPA3-expressing 
sf9  cells  [1].  These  same  modifications  resulted  in  inactive 
compounds  when  the  LPA3  splice  variant,  HOFNH30,  was 
expressed  in  RBL-2H3  cells  [3].  In  whole  animal  assays  in 
rats,  replacement  of  X  with  choline,  ethanol,  propanol,  or 
methanol  results  in  hypotension,  rather  than  the  hypertension 
normally  induced  by  LPA  [13].  Several  combined  modifica¬ 
tions  in  the  X  and  Q  regions,  including  the  palmitoyl 
dioxazaphosphocanes  DOXP-H  and  DOXP-OH  (Fig.  1), 
were  tested  for  their  effects  on  platelet  aggregation  [14]. 
Interestingly,  DOXP-OH,  with  its  ionizable  phosphoiyl 
group,  was  the  only  structure  to  show  LPA-Iike  effects.  On 
the  other  hand,  DOXP-H  inhibited  LPA-stimulated  platelet 
aggregation.  Platelets  have  been  shown  to  express  all  three  of 
the  LPA-selective  members  of  the  EDG  family  (LPAi/EDG2, 
LPA2/EDG4  and  LPA3/EDG7)  [15].  Thus,  it  is  still  unclear  if 
DOXP-H  might  be  a  general  LPA  receptor  antagonist,  or 


show  some  specificity  toward  one  or  two  LPA  receptors.  One 
modification  that  is  well  tolerated  by  all  three  receptors  is 
cyclization  of  the  phosphate  with  the  2-hydroxyl  of  natural 
LPA  to  form  cyclic  phosphatidic  acid  (cPA,  Fig.  1)  [1,12]. 
This  modification  maintains  both  the  approximate  size  and 
the  charge  of  the  headgroup.  These  experimental  studies,  as 
well  as  the  recent  experimentally  validated  theoretical  model 
of  the  LPA]  complex  with  LPA  [16]  all  point  to  a  critical 
interaction  between  the  anionic  phosphate  group  of  LPA  and 
its  receptors. 

The  most  extensively  examined  series  of  LPA  modifica¬ 
tions  have  been  in  the  region  labeled  R1  in  Fig.  1.  Assays  of 
Ca2+-mobilization  in  LPA3-expressing  insect  sf9  cells  [1,17] 
and  GTP-y-35S  binding  in  LPA3-expressing  HEK293T  cells 
[2]  both  indicate  that  unsaturated  acyl  chains  are  strongly 
preferred  over  saturated  acyl  chains,  whereas  LPA!  and  LPA2 
expressed  in  the  same  cell  types  do  not  show  such  selectivity. 
This  strong  preference  of  LPA3  for  unsaturated  acyl  chains  as 
well  as  the  expression  of  LPA3  in  the  heart  [1 ,2]  indicates  that 
the  cardiovascular  effects  in  rats  that  are  more  strongly 
induced  by  unsaturated  acyl  chains  [13]  are  likely  mediated 
by  the  LPA3  receptor.  All  three  receptors  show  considerably 
decreased  or  no  activity  when  the  acyl  chain  length  is  reduced 
below  16  [2,17].  A  combined  modification  ofR1  and  R2  in  the 
form  of  2-acyl  LPA  has  also  been  tested  against  LPA2  and 
LPA3-expressing  sf9  cells.  LPA3  shows  a  strong  preference 
for  2-acyl  LPA  species  over  the  corresponding  1-acyl  LPA 
species,  while  LPA2  shows  no  such  preference  [17]. 

Limited  modifications  to  the  region  labeled  Q  in  Fig.  1 
have  been  examined  in  cells  overexpressing  single  LPA 
receptors.  In  particular,  studies  with  HEK293T  cells  ex¬ 
pressing  either  LPAi,  LPA2  or  LPA3  demonstrate  that  the 
LPA]  and  LPA2  receptors  both  accept  ethanolamine  as  a 
substitute  for  glycerol,  in  fact,  both  are  more  responsive  to 
the  ethanolamine  derivatives  [2].  LPA3,  on  the  other  hand, 
very  poorly  responds  to  the  ethanolamine  compounds. 


3.  LPA  receptor  antagonists 

Within  the  last  year,  two  separate  groups  have  reported  the 
characterization  of  selective  antagonists  for  LPA  receptors. 
Fischer  et  al.  reported  the  testing  of  a  selective  and  compet¬ 
itive  inhibitor  of  LPA3,  diacylglycerol  pyrophosphate 
(DGPP)  8:0  (Fig.  1).  DGPP  8:0  inhibits  LPA-induced  Ca2+ 
responses  in  RH7777  cells  expressing  LPAi  and  LPA3,  with 
Ki  values  of  6600+680  and  106+28  nM,  respectively  [18]. 
Fig.  2  demonstrates  this  selective  inhibition  of  the  LPA- 
induced  Ca2+  responses  in  stably  transfected  RH7777  cells 
expressing  the  LPA3  and  LPAb  but  not  in  either  stably  or 
transiently  transfected  RH7777  cells  expressing  the  LPA2 
receptor.  Heise  et  al.  reported  that  an  N-acyl  ethanolamide 
phosphate  (NAEPA)  with  a  benzyl-4-oxybenzyl  group  at  the 
2-position  of  the  ethanol  backbone  (Fig.  1)  selectively 
inhibited  LPA-stimulated  GTP[y35S]  binding  in  stably  trans¬ 
fected  LPA]  and  LPA3-expressing  RH7777  cells,  with  Kt 
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LPAi  -  stable  LPA2  -  stable 


LPA  +  DGPP 


Fig.  2.  DGPP  8:0  inhibits  the  Ca2+ responses  elicited  by  LPA  in  LPA]  andLPA3,  but  not  LPA2  expressing  RH7777  cells.  RH7777  cells,  expressing  LPAb  LPA2,  or 
LPA3  receptors,  were  exposed  to  a  mixture  of  100  nM  (Panels  A,  C  and  D)  or  1  pM  (Panel  B)  LPA  18:1  and  1  pM  (Panels  A,  C  and  D)  or  10  pM  (Panel  B)  DGPP 
8:0.  Control  cells  were  exposed  to  1 00  nM  (Panels  A,  C  and  D)  or  1  pM  (Panel  B)  LPA  1 8 : 1 .  Representative  Ca2+  responses  are  shown  for  stable  LPAi  (Panel  A), 
LPA2  (Panel  B),  and  LPA3  (Panel  C)  expressing  cells,  or  cells  transiently  expressing  LPA2  (Panel  D).  Stable  LPA2  cells  had  a  higher  EC50  for  LPA  as  compared  to 
cells  transiently’  expressing  the  same  construct  Nevertheless,  the  high  10  pM  concentration  of  DGPP  8:0  was  still  ineffective  in  inhibiting  the  response. 


values  of  137  and  428  nM,  respectively  [19].  Since  the  assays 
used  to  characterize  these  antagonists  were  not  the  same,  the 
Ki  values  cannot  be  compared  to  determine  the  relative 
efficacy  of  the  antagonists.  It  is  interesting  to  note,  however, 
that  neither  antagonist  shows  activity  against  LPA2/EDG4, 
and  that  they  show  opposite  preferences  for  LPAi  and  LPA3. 

These  antagonist  structures  share  a  significant  feature  with 
the  known  SAR  for  the  LPA  receptors.  The  polar  headgroup 
of  both  structures  is  anionic,  either  phosphate  or  pyrophos¬ 
phate.  However,  both  of  these  structures  contain  two  hydro- 
phobic  chains,  a  considerable  change  from  the  typical  agonist 
at  these  receptors.  One  hydrophobic  chain  in  each  structure  is 
relatively  short,  8  carbons  in  the  case  of  DGPP  and  approx¬ 
imately  10  carbons  long  in  the  case  of  the  benzyl-4-oxy- 
benzyl  group.  The  structures  differ  considerably  in  the  com¬ 
position  of  the  second  lipid  chain,  which  includes  8  carbons 
in  the  case  of  DGPP  8:0  but  18  carbons  in  the  case  of  the 
NAEPA  analog.  The  presence  of  one  short  hydrophobic  chain 
seems  to  be  important  as  neither  DGPP  1 8:0  norphosphatidic 
acid  18:0  were  antagonists  at  the  LPA  receptors  [18]. 

4,  Receptor  structure  and  amino  acids  critical  for  ligand 
recognition 

The  intriguing  similarities  and  differences  between  the 
known  agonists  and  antagonists  of  the  LPA  receptors  led  us 
to  apply  molecular  modeling  techniques  to  explore  the  basis 
for  both  the  observed  antagonist  selectivity  for  LPAi  and 
LPA3,  and  for  agonism  versus  antagonism. 


4.1.  LPA  receptor  modeling 

Models  of  the  LPA  receptors  were  developed  by  homol¬ 
ogy  to  our  previously  validated  SI  Pi  model  [20].  The 
validated  SIP!  model  is  consistent  with  experimental  studies 
exploring  the  impact  of  site-directed  mutations  on  receptor 
activation,  and  thus  these  models,  referred  to  as  the  ‘A* 
models,  represent  the  active  conformation  of  these  receptors. 
Since  antagonists  should  have  equal  affinity  for  the  inactive 
conformation,  we  also  developed  models  of  each  receptor 
based  on  the  recently  published  bovine  rhodopsin  crystal 
structure,  which  was  characterized  in  its  dark,  inactive  state, 
termed  the  T  models  [21].  These  T  models,  since  they  are 
based  on  the  only  experimental  structure  of  a  G  protein- 
coupled  receptor,  are  likely  to  accurately  represent  regions 
outside  the  transmembrane  domain.  This  is  not  the  case  for 
models  based  on  the  SlPj  model,  which  was  based  on  a 
template  structure  that  did  not  have  atomic  coordinates  out¬ 
side  the  transmembrane  domain.  The  validation  of  the  SI  Pi 
model  additionally  includes  only  critical  amino  acid  contacts 
with  S  IP,  all  of  which  are  within  the  transmembrane  domains. 
All  receptor  models  were  developed  and  optimized  according 
to  the  procedures  we  have  applied  to  the  development  of  our 
validated  SlPi/EDGl  and  LPA:/EDG2  models  [16,20]. 

Ballesteros  and  Weinstein  [22]  developed  a  unique 
nomenclature  for  identifying  residues  that  are  in  correspond¬ 
ing  locations  in  different  GPCRs.  This  nomenclature  uses 
the  most  conserved  residue  of  a  helix  and  identifies  it  with 
the  helix  number  and  then  the  extension  of  50  (example: 
W4.50  is  the  most  conserved  residue  in  the  fourth  helix). 
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Fig.  3  identifies  each  helical  reference  position  on  an 
alignment  of  bovine  rhodopsin  and  the  LPA  receptor  sub¬ 
family.  All  residues  in  the  transmembrane  domain  can  be 
numbered  relative  to  the  helical  reference  position.  This 
nomenclature  is  particularly  useful  for  identifying  amino 
acids  that  are  in  the  same  position  within  the  transmembrane 
domains  of  different  receptors  and  so  will  be  used  through¬ 
out  this  text  in  order  to  better  compare  the  LPA1?  LPA2  and 
LPA3  complexes. 

4.2.  Docking  studies 

Docking  computations  were  performed  to  determine 
interactions  between  ligand  and  receptor.  The  resulting 
complexes  were  then  used  to  determine  whether  differences 


exist  between  the  members  of  the  LPA  receptor  family  that 
explain  die  experimentally  observed  differences  in  affinity 
for  DGPP.  Two  different  docking  methods  were  utilized  to 
fully  explore  die  interactions  between  DGPP  and  the  LPA 
receptors  and  to  quantitatively  correlate  computed  docking 
energies  with  the  experimentally  observed  trend  in  selectiv¬ 
ity.  The  Tabu  algorithm  [23]  in  the  MOE  computer  program 
(version  2001.01,  Chemical  Computing  Group,  Montreal) 
with  random  ligand  starting  positions,  was  first  used  to 
generate  up  to  50  complexes  for  the  appropriate  receptor 
models  and  either  DGPP  or  1 8: 1  LPA.  The  resulting  putative 
complexes  were  evaluated  as  in  our  previous  studies  for  the 
presence  of  electrostatic  interactions  between  the  receptor 
and  ligand  [20,24,25].  Autodock  3.0  [26]  was  additionally 
used  to  generate  100  complexes  of  DGPP  with  the  ‘I*  models 
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s  s  8  as 

^  oi  cvi  on 

51  QFPI HFLTLYVTVQHKKLRTPLNY I LLHLAVADLFMVFGGFTTTLYTSLHGYFVFGPTOCNLEOF 
63  I MLANLLVMVA I YVNRRFHFP I YYLMANLAAADFFA6LAYFYLMFNTGPNT  ■(§RLTVSTWLl|SQG 
46  VLLTMLLVIAAIASNRRFHQPIYYLLGNLAAADLFAGVAYLFLMf0TGPRT.ARLSLEGKVFlRaG 
44  | FFSNSLV I AAV I KNRKFHFPFYYLLANLAAADFFAG I AYVFLMFNTGPVS  ■®TLTVNRWFljgOG 


s  f 

116  FATLGGE I ALWBLVVLA I ERYVVVCKPMSH FRF6EHHA I MGVAFTWVMALACAAPP  LVQWSRY I P 
127  LI DT8LTA8VANLLA I A I ERH I TVFRMQLHTRMSNRRVVVV I W IWTMA I VMGA I PSVGWNC I C • 
110  LLDTSLTASVATLLAI AVERHRSVMAVQLHSRLPRGRWML I VGVWVAALGLGLLPAHSWHCLC- 
108  LLDSSLTASLTNLLV I AVERHMS I MRMRVHSNLTKKRVTLL I LLVWA I A I FMGAVPT  LGWNCLC • 


S? 

181  FGMQCaCQ I DYYTPHEETNNE3FV I YMFWHF II  PL I V I FFCYGQLVFTVKEAAA • • -QQQESAT 


191  . Dl  ENCSNMAPLYSDSYLVFWAI  FNLVTFWMWLYAHI  FG •  YVRQRTMRMSRHSSGPRR 

174 . ALDRC8RMAP  L  LSRSYL AVWALSS  L  LV  F  L  LMVAVYTR I FF ■ YVRRRVQRMAEHVSCHPR 

172 . Hi  SACSSLAP I  YSRSYLVFWTVSNLMAFL  IMVVVYLRI  YV  ■  YVKRKTMVL8PHTSGS  I S 


8  8  S 

243  TOKAEKEVTRMV I  I MV I AFL I CWLPYAGVAFY I FTHQGSDFGP I FMT I PAFFAKTSAVYNPV I Y I 
249  NRDTMHS LLKTVV I VLGAF I  1 CWTPGLVLLLLD. VCCPQCDVLAY^KjFFLLLAEFNSAMNP I IYS 
232  YRETTLSLVKTVVI I LGAFVVCWTPGQWLLLDGLGCESCNVLAVEKYFLLLAEANSLVNAAVYS 
230  RRRTPMKLMKTVHUTVLGAFVVCWTPGLWLLLDGLNCRQCGVQHVI^gliyFLLLALLNSVVNP I  IYS 


Rhod  308  IWNKQFRNCIIIVTTLCCGK . NPLGDDEASTTVSKTETSQ*  >VAPA 

LPA,  313  YRDKEMSATFRQI LCCQR- • SENPTGP ■ ■ ■ TEGSDRSASSL . NHT I LAGVHSNDHSW 

LPAa  297  CRDAEMRRTFRRLLCCAC* • LRQSTRESVHYTSSAQGGASTRIMLPENGHPLMTPPFSYLELQRY 
LPAs  295  YKDEDMYGTMKKMICCFSQENPERRPSRI PSTVLSRSOTG* - SQY I  EDS  I SQGAVCNKSTS 


LPAa  360  AASNKSTAPDDLWVLLAQPNQ 


Fig.  3.  Aligned  sequences  for  the  LPA  family  and  thodopsin.  From  top,  bovine  rhodopsin,  human  LPA,,  human  LPA2  and  human  LPA3.  Helical  reference 
positions  are  indicated  above  the  appropriate  position,  as  are  specific  positions  mentioned  in  the  text  Transmembrane  helices  including  extensions  beyond  the 
membrane  observed  in  the  bovine  rhodopsin  crystal  structure  are  underlined.  Boxes  highlight  residues  that  docking  studies  indicate  are  important  in  binding  file 
pyrophosphate  group  of  DGPP  8:0. 
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of  each  LPA  receptor  in  order  to  more  completely  explore  die 
binding  site  found  by  the  Tabu  algorithm  in  the  vicinity  of  the 
extracellular  loops,  as  well  as  to  provide  a  more  quantitative 
reflection  of  the  relative  affinity  of  DGPP  for  each  receptor. 


The  best  complex  for  each  receptor  was  then  energy  mini¬ 
mized  with  the  MMFF94  forcefield  to  a  gradient  of  0.0 1  kcal/ 
mol  A  to  allow  the  receptor  to  conform  to  the  presence  of 
DGPP.  DGPP  was  then  removed  and  100  additional  com- 


Fig.  4.  Models  of  LPA  receptors  with  18:1  LPA.  Proteins  are  shown  as  ribbon  models  with  18:1  LPA  as  a  space-filling  model.  Residues  forming  key 
interactions  with  the  polar  headgroup  of  18:1  LPA  are  shown  as  stick  models.  The  transmembrane  domains  are  numbered  in  Panel  A  using  Roman  numerals, 
the  orientation  of  these  domains  in  Panels  B  and  C  is  the  same.  LPAj  (Panel  A),  LPA2  (Panel  B)  and  LPA3  (Panel  C)  complexes  with  18:1  LPA  viewed  from 
within  the  cell  membrane.  Close-up  of  interactions  between  LPAj  (Panel  D),  LPA2  (Panel  E)  and  LPA3  (Panel  F)  and  the  polar  headgroup  of  18:1  LPA  viewed 
from  the  same  perspective  as  in  Panels  A,  B  and  C.  (Panel  G)  Close-up  of  interactions  between  LPA3  and  the  polar  headgroup  of  1 8:1  LPA  viewed  from  the 
extracellular  surface  (rotated  approximately  90°  from  the  orientation  shown  in  Panel  F). 
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plexes  were  generated  with  Autodock  3.0  in  order  to  reeval¬ 
uate  the  binding  energy  after  allowing  the  receptor  model  to 
adapt  to  the  presence  of  the  ligand.  If  the  position  found  prior 
to  minimization  was  not  found  in  the  100  additional  com¬ 
plexes  generated,  the  docking  run  was  repeated.  Once  a 
docking  run  succeeded  in  finding  the  pre-minimized  posi¬ 
tion,  the  complex  with  the  lowest  docked  energy  (whether 
the  ligand  was  in  the  pre-minimized  position  or  not)  was 
selected  as  the  best  complex. 

4.3.  18:1  LPA  complexes 

Using  18:1  LPA  and  the  ‘A’  models  of  LPA  b  LPA2  and 
LPA3,  the  complexes  generated  during  the  docking  studies 


were  evaluated  to  determine  the  amino  acids  involved  in 
agonist  binding.  On  examination  of  the  LPA  receptor  family, 
a  clear  and  consistent  picture  of  LPA  binding  emerged.  Three 
residues  are  predicted  to  be  necessary  for  binding  LPA: 
R3.28,  Q3.29  and  K/R7.36,  Fig.  4.  An  additional  residue, 
K7.35,  may  be  important  for  LPA  binding  to  LPA3.  In  LPAb 
mutation  of  Q3.29  (Q125)  to  A  completely  abolished  binding 
of  LPA  [16].  This  prediction  is  also  consistent  with  the  SIP 
docking  studies  with  SlPj  [20,24,25].  These  studies  exper¬ 
imentally  validate  the  importance  in  SIPi  of  R3.28,  E3.29 
and  R7.34  for  binding  of  S  IP.  In  the  alignment  of  the  LPA  and 
SIP  receptor  families,  positions  3.28  and  3.29  are  conserved 
across  the  two  families,  but  die  position  of  the  third  residue  is 
not  as  strictly  conserved  in  the  SIP  receptor  family. 


Fig.  5.  Models  of  LPA  receptors  with  8:0  DGPP.  Proteins  aTe  shown  as  ribbon  models  with  8:0  DGPP  as  a  ball  and  stick  model.  Residues  forming  key 
interactions  with  the  polar  headgroup  of  8:0  DGPP  are  shown  as  stick  models  in  Panels  D-F.  The  transmembrane  domains  are  numbered  using  Roman 
numerals  in  Panel  A,  the  orientation  of  these  domains  is  the  same  in  Panels  B  and  C.  LPAi  (Panel  A),  LPA2  (Panel  B)  and  LPA3  (Panel  C)  complexes  with  8:0 
DGPP  viewed  from  within  the  cell  membrane.  Close-up  of  interactions  between  LPA!  (Panel  D),  LPA2  (Panel  E)  and  LPA3  (Panel  F)  and  the  polar  headgroup 
of  8:0  DGPP  viewed  from  the  same  perspective  as  in  Panels  A  B  and  C. 
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4.4 .  DGPP  8:0  complexes 

The  complexes  generated  by  the  MOE  program  were 
evaluated  for  interactions  between  die  receptor  and  the 
ligand  based  on  the  number  of  hydrogen  bonds.  The 
complex  with  the  most  hydrogen  bonds  was  used  to 
determine  the  optimal  region  of  the  receptor  structure  to 
more  thoroughly  explore  with  the  Autodock  program. 
Docked  energies  from  the  Autodock  program  were  utilized 
to  select  the  best  complex  in  order  to  compare  interactions 
across  the  LPA  receptor  family. 

Evaluation  of  the  complexes  between  DGPP  and  the  T 
model  of  the  DGPP-selective  LPA3  receptor  produced  multi¬ 
ple  entries  that  occupied  the  same  volume  at  the  extracel¬ 
lular  end  of  the  transmembrane  helical  bundle,  yet  within 
the  extracellular  loops.  The  complex  with  the  lowest  docked 
energy  (-12.4  kcal/mol)  is  shown  in  Fig.  5.  The  polar  head 
group  forms  hydrogen  bonds  with  residues  K95  (extracel¬ 
lular  loop  1),  R105  (3.28),  N172  (extracellular  loop  2)  and 
R276  (7.36).  The  complex  of  LPA  and  LPA3  indicates  that 
R3.28  and  R7.36  are  necessary  for  agonist  binding  and  their 
involvement  in  binding  both  the  agonist  and  antagonist  is 
consistent  with  the  experimental  result  that  DGPP  is  a 
competitive  antagonist  of  LPA  action  at  LPA3.  The  hydro¬ 


phobic  tails  of  DGPP  interact  entirely  with  residues  at  the 
extracellular  end  of  the  transmembrane  domains  and  with 
residues  in  the  extracellular  loops.  One  hydrophobic  chain  is 
oriented  toward  transmembrane  helices  IV  and  V.  The 
second  hydrophobic  chain  is  oriented  toward  transmem¬ 
brane  helix  I  (Fig.  6). 

Complexes  of  the  LPAi  receptor,  model  T,  with  DGPP 
showed  geometric  features  similar  to  the  best  complex  of  the 
higher  affinity  LPA3  receptor  complex,  but  with  a  poorer 
docked  energy  (—11.9  kcal/mol).  Fig.  5  shows  the  best 
complex  between  DGPP  and  LPAi.  The  polar  headgroup  of 
DGPP  forms  hydrogen  bonds  with  R36  (amino-terminus), 
R114  (extracellular  loop  1),  R124  (3.28)  and  K294  (7.36). 
The  complex  of  LPA  and  LPAi  predicts  that  R3.28  and 
K7.36  are  necessary  for  agonist  binding  and  their  interaction 
with  the  antagonist  is  consistent  with  the  experimental  result 
that  DGPP  competitively  antagonizes  LPA  action  at  LPAi. 
The  hydrophobic  tails  of  DGPP  interact  entirely  with 
residues  at  the  extracellular  end  of  the  transmembrane 
domains  and  with  residues  in  the  extracellular  loops.  The 
hydrophobic  tails  in  this  complex  are  both  oriented  in  the 
direction  of  helices  IV- VI  (Fig.  6). 

The  best  DGPP  complex  generated  with  the  DGPP- 
unresponsive  LPA2  has  a  significantly  higher  docked  energy 


Fig.  6.  Overlay  of  DGPP  complexes  with  LPAl5  LPA2  and  LPA3.  For  clarity,  ribbons  are  shown  only  for  the  extracellular  loops  and  the  first  helical  turn  into  the 
membrane  for  LPA2,  ribbons  for  the  other  two  receptors  are  analogous.  The  transmembrane  helical  domains  are  numbered  with  Roman  numerals.  DGPP 
structures  are  shown  as  stick  models  with  the  phosphorous  atoms  as  larger  balls.  The  green  DGPP  structure  is  the  optimal  position  in  LPAlt  blue  shows  the 
optimal  position  in  LPA2  and  orange  shows  the  optimal  position  in  LPA3. 
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(—9.4  kcal/mol)  than  the  DGPP-selective  LPA3  complex 
(-12.4  kcal/mol)  or  the  moderately  responsive  LPAi  com¬ 
plex  (-11.9  kcal/mol).  These  predicted  differences  agree 
with  the  experimental  observation  that  LPA  action  at  LPA2  is 
not  inhibited  by  DGPP  while  DGPP  demonstrates  a  106  nM 
Kt  at  LPA3  and  a  6600  nM  Kj  at  LPA!  [18].  Only  three 
residues  were  found  to  form  hydrogen  bonds  with  DGPP, 
H91  (2.63),  R107  (3.28)  and  K278  (7.36)  (Fig.  5).  Two  of 
these,  R3.28  and  K7.36  are  predicted  to  interact  with  LPA. 
The  hydrophobic  chains  of  DGPP  in  this  complex  are 
oriented  toward  transmembrane  helices  I  and  VII  (Fig.  6). 


5.  Implications  of  selective  LPA  receptor  antagonists 

The  LPA i -selective  antagonist,  benzyl-4-oxybenzyl 
NAEPA,  and  the  LPA3-selective  antagonist,  DGPP  8:0,  are 
invaluable  pharmacological  tools  for  researchers  examining 
the  signaling  pathways  of  the  LPA  receptors  in  numerous  cell 
types.  The  endogenous  co-expression  of  multiple  LPA- 
responsive  receptors  in  a  wide  variety  of  cells  has  previously 
been  a  significant  challenge  to  researchers  wishing  to  exam¬ 
ine  the  biology  of  LPA  without  the  complication  of  potential 
artifacts  due  to  receptor  overexpression. 

The  modeling  results  presented  herein  indicate  that  the 
molecular  basis  for  the  observed  competitive  antagonism  by 
these  compounds  is  due  to  a  common  interaction  site  for  the 
anionic  headgroup  of  both  agonists  and  antagonists.  The 
anionic  phosphate  group  of  both  DGPP  and  LPA  is  anch¬ 
ored  to  the  receptor  through  interactions  with  a  common  set 
of  cationic  amino  acids  in  transmembrane  helices  III  and 
VII.  However,  a  very  different  interaction  exists  between 
DGPP  and  the  extracellular  loops.  Preliminary  docking 
studies  with  benzyl-4-oxybenzyl  NAEPA  indicate  that  this 
antagonist  interacts  with  the  extracellular  loops  in  a  manner 
similar  to  DGPP.  These  results  suggest  that  antagonists  can 
be  developed  that  completely  diverge  from  the  phospholipid 
structures  explored  thus  far  against  these  receptors. 

Finally,  LPA  receptor  antagonists  can  now  be  explored  as 
potential  therapeutic  lead  compounds.  The  importance  of 
LPA  in  angiogenesis  [27,28],  inflammation  [29],  cancer 
invasiveness  [30]  and  atherosclerosis  [31,32]  demonstrates 
but  a  few  of  the  potential  processes  that  might  be  targeted  by 
LPA  receptor  agonists  and  antagonists. 
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|  MOLECULAR  MODELING  OF  LYSOPHOSPHATIDIC  ACID  RECEPTOR 
i  ANTAGONISTS.  Vineet  M.  Sardar1,  Tamas  Virag2,  David  J.  Fischer2,  Don 
j  Elrod3,  Debra  L.  Bautista1,  De-an  Wang2,  Nora  Nusser2,  KazUaki  Yokoyama2, 

;  Daniel  L.  Baker2,  Duane  D.  Miller3,  Gabor  Tigyi2,  and  Abby  L.  Parrill1.  (1) 
Department  of  Chemistry  and  Computational  Research  on  Materials  Institute, 
j  University  of  Memphis,  3744  Walker  Suite  /  213,  Memphis,  TN  38152-6060, 
Fax:  901-678-3447,  vsardar@memphis.edu,  (2)  Department  of  Physiology  and 
1  Biophysics,  University  of  Tennessee,  (3)  Department  of  Pharmaceutical 
j  Sciences,  University  of  Tennessee 

Lysophosphatidic  acid  (LPA)  elicits  a  variety  of  responses  including  mitogenesis, 
|  cytoskeletal  changes,  activation  of  Ca2+  transients,  and  effects  on  apoptosis. 

I  These  responses  are  elicited  via  LPA1/EDG2,  LPA2/EDG4,  and  LPA3/EDG7  G 
protein-coupled  receptors.  These  receptors  are  members  of  the  endothelial 
differentiation  gene  family.  In  order  to  understand  the  physiological  significance 
|  of  LPA  highly  selective  antagonists  are  necessary.  Recently,  dioctyl  glycerol 
pyrophosphate  (DGPP)  and  fatty  alkyl  phosphate  (FAP)  were  shown  to  be  potent 
and  selective  antagonists  towards  LPA3  receptor.  We  have  docked  DGPP  and 
FAP  in  our  LPA,,  LPA2,  and  LPA3  receptor  models  and  our  docked  energies 
agree  with  the  observed  trend  in  inhibition  constants  (K;).  The  docked  positions 
of  the  antagonist  relative  to  the  agonist  overlap  in  the  position  of  the  polar  head 
group,  but  diverge  in  the  favored  position  of  the  hydrophobic  tail(s). 


